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Abstract 

An effective action approach to Kohn-Sham density functional theory is used to illustrate how 
the exact Green's function can be calculated in terms of the Kohn-Sham Green's function. An 
example based on Skyrme energy functionals shows that single-particle Kohn-Sham spectra can 
be improved by adding sources used to construct the energy functional. 
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The Skyrme-Hartree-Fock approach to nuclear properties has had wide success in repro- 
ducing bulk properties of nuclei across the periodic table [H 0, 0, 0, IE E( • The interpretation 
of the Skyrme formalism as an approximate implementation of Kohn-Sham density func- 
tional theory (DFT) |7] implies that certain observables (energy per particle, densities) can 
be calculated reliably, but these do not include single-particle quantities. Only for the bulk 
observables can the DFT framework accommodate all correlations in principle (if not in 
practice because of the limited form of the energy functionals actually used) || 0, Eof . 
Nevertheless, single-particle energies and wave functions from Skyrme and other DFT-like 
formalisms are also regularly used. 

In this letter, we illustrate how to extend the effective action approach to Kohn-Sham 

DFT H E3 

Il4j to calculate the full single-particle Green's function in terms of Kohn- 



Sham Green's functions at the same level of approximation. Our discussion directly adapts 
the extension described in the context of Coulomb systems in Refs. HI El. This connection 



between Green's functions helps to clarify both some misconceptions and limitations of the 
Kohn-Sham approach, and suggests how to improve calculations of single-particle properties. 
At first, we consider functionals of the fermion density only, and then compare to generalized 
functionals that also depend on the kinetic energy density to illustrate the effect of additional 
sources. 

We introduce a generating functional in the path integral formulation with a Lagrangian 
C supplemented by a local c-number source J(x) coupled to the composite density operator 
as m Ref. [[1, but add a non-local c-number source £(x,x') a /3 coupled to tp a (x)ip^(x') , 



Z[J £l — e iW{J,€] _ / JJ^jjJJ^ e i I d ' lx \ C + J{ x )i'l t { x )^c l (x) + ^d^x 1 ■tp a (x)i(x,x') a , 3 -tp^(x')] /j\ 



where a and (3 are spin indices and summation o f rep eated indices is implied. (We generalize 
below to an additional local source, as in Ref. [13].) For simplicity, normalization factors 
are considered to be implicit in the functional integration measure 17, 0]. As a specific 
example, we will use the effective field theory (EFT) Lagrangian appropriate for a dilute 
Fermi system 1^|, but the discussion can be adapted to any system for which a hierarchy 



of approximations can be defined. 

The fermion density in the presence of the sources J and £ is 

Note that the sources here are time dependent, in contrast to the more limited discussion 
with static sources in Ref. 0]; however, the generalization of the formalism is direct. A 
functional Legendre transformation from J to p, which takes us from W to the effective 
action T, produces an energy functional of the density, which is minimized at the exact 
ground-state density for time-independent sources. 1 The inversion method 0, 0] carries 
out this inversion order-by-order in a specified expansion; an EFT expansion was used in 
Refs. 3 and [l3j]. At the end, one sets J(x) and £(x, x') to zero. (Although we are unaware 



of any general problems, we have not excluded the possibility of complications in making 
the inversions with time-dependent sources.) 



Note that the energy functional is only obtained once £ is set to zero. 
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Solving the zeroth-order system defines the Green's function G ks (x, x') a p of the Kohn- 
Sham non-interacting system in the presence of £(x,x') a p, the Kohn-Sham potential Jq(x), 
and an external potential v{x). This Green's function satisfies 

J d 4 z [G ks (x, z)]^Gte(z, x%p = 5 a p5 4 (x - x) , (3) 



or 

d 4 z (id t + ^ - v(x) + J {x)j 5 aJ S 4 (x -z)- £(z, x) ja G ks (z, x') 1 p = 5 a p5 4 (x - x') , 

(4) 

with appropriate finite-density boundary conditions (one could also introduce a chemical 
potential). Note that G ks doesn't take a simple form in terms of orbitals [see G° s in Eq. 
until we set £ = and restrict ourselves to time independent J . 

Functional derivatives of W with respect to £(x,x') gives the two-point function in the 
presence of the sources, 

The exact ground-state Green's function is obtained by setting £ = J = after taking 
derivatives. The key results we will need to evaluate Eq. (JBJ) in terms of Kohn-Sham quanti- 
ties were given in Refs. [TBI Il6| (we follow their notation for the most part) and are rederived 
here. First, functional derivatives with respect to £ of W and F are directly related, where 

r[p,£] = W[J,Z]- Jd 4 y.J(y)p(y) (6) 

is the effective action. Namely, the functional derivative with respect to £ of this equation 
yields (spin indices are suppressed) 

p = \6&$) // A \~sw) t [wBv ~S dAy Gec^O ) P p{y) ' 

from which the last two terms cancel, leaving 

6W[J,£] \ ( $?\p,i\ 



8i{x } x') a pJ j \8£{x,x') a p 
(Here and below we repeatedly apply the functional relations 



(7) 
(8) 



\6U P \8tJj \5JJs\SU p \8pJt\5t 

where F = F[J,£] and arguments and integrals are implied.) Equation (|SJ) is a special case 
of a general result for Legendre transformations proved in Ref. |20J. 

Next, this relation applied to the zeroth-order (Kohn-Sham) system yields the Kohn- 
Sham Green's function, 

Ura), = im^u) Ja = ,G ^ xU ■ (10) 
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We divide the full effective action into zeroth-order and interacting pieces, 

r[M] = r [p,£] + r int [ P ,£] . 

Since r int [p, £] depends on £ only through Gks, 



5T int [p,£] / SG ks (y,y' 



5G ks (y,y') 5l \ 5£(x,x') a p 



)6 n 



dydy' . 



(11) 



(12) 



The second half of the integrand can be rewritten 



8G ks (y,y') Sl \ ( SG ks (y,y') Sl 



5£(x,x r 



aj3 



8^{x,x') aP J Jo 
$G ks (y,y') s ~f 



+ 



5G ks (y,y') s ^ \ ( SJ (z) 
8G kB {y : y') Sl 



5£(x,x') a p 



dz 



G ks {x,y') ai G ks (y,x')si3 

$G ks (y,y') S ry 



Sp(z) 



dz 



Jo 



+ 1 



8p(z) 



G ks (x,z) aX G ks (z,x')x/3dz . 



(13) 



The second line follows by applying Eq. with F — > G ks and simplifying. The functional 
derivatives in the second line can be evaluated by using the expression for G ks in terms of 
the noninteracting generating functionals. Thus, 
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SG ks (y,y') Sl 

5£{x,x') a p J j S^(x,x') 



la/3 



1 6Z [J ,£] 
Zo5£(y,y') 



1 6Z [J ,£] 
Z 8£(x,x') a p 



1 6Z [J ,£] 
Z 5^(y,y')sj 



1 S 2 Z [J ^] 



Z 5£,(x,x') a p 5£(y,y') s ~f 

= i*(T[Mx)1>iW)]) (TMvWW) -i 2 <T[M*HW)MvWW) 
= (-if(T[MM(y')}) (T[Mv)1>lW]) 

= G ks (x,y') ay G ks (y,x')si3 , (14) 

where we've applied Wick's theorem to the noninteracting system to go from the third line 
to the fourth line, and 



Sp(z) 



5£,(x,x') a/3 J Jo V ^(x,x') aP ) Jo 



6G ks {z,z + )ss 



-iG ks (x, z) aS G ks (z,x')s/3 ■ (15) 



Alternatively, we can expand 5(G ks G k ^) / 6£ = and use SG k ^ /<5£ = —1. 
Substituting Eq. (JT3j) back into Eq. (JT2J), we find that 



5£(x,x') aP 



G ks (x, y') ai ^tt7^4- G ks(y, x') Sl 3 dy dy' 



+ i / G ks (x,y) aX 



SG ks (y,y') Sy 

5r int [p, c] 



8p(y) 



G ks (y,x') X f3 dy 



(16) 
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FIG. 1: Equation for the full propagator in terms of the Kohn-Sham Green's functions and self- 
energy. 



Equations (fTUJl and (jltij). together with 



iG(x,x') a p 



5£{xix') a p) j \8^(x,x') a p / „ 
imply a Dyson equation for the exact Green's function: 



8^(x,x') aP/ p 



8T int [p,£] 

8£(x,x') a p 



, (17) 



G(x, x') a p = G ks (x, x') a p + / G ks (x, y')aj ^ks(y', y) 7 s G ks (y, x')sp dy dy' 



(18) 



which defines a self-energy S ks as 



£ks(2/',y) 



7# 



i 5r int [p,e] 



+ 



flr int [p, f] 



i8G kB {y,y') Sl V M?/) 



S-yS 8 A (y' - y) 



(19) 



In the second line, the self-consistent Kohn-Sham potential Jo is equal to 8Y int /8p only when 
we set J = 0| . Neither E ks nor S ks is the conventional self-energy, which is built from 
non- interacting (rather than Kohn-Sham) Green's functions. We can obtain E ks (y',y) at 
the diagrammatic level by opening each G ks line in turn in a given Feynman diagram for 
r int . It consists of the same diagrams as the conventional one-particle-reducible self-energy, 
but with the fermion lines given by Gks rather than the non-interacting Green's function 
(which includes only the external potential). 

Now consider applying these equations with £ = J = after taking functional derivatives; 
we denote the Kohn-Sham Green's function in this case as G° s . For simplicity we will 
consider spin-independent interactions, so that the Green's functions and self-energies are 
diagonal in spin. Kohn-Sham orbitals arise as solutions to 



2M 



v x 



J (x)] ^ fc (x) = £fc^fc(x) , 



(20) 



where the index k represents all quantum numbers except for the spin |14j. The decompo- 
sition of G ks (x, x') a p = SapG^x, x') in terms of these orbitals is 14] 

iG^KxY) = J> fc (x) ^(x') e-**W[e(t - 6{e k - e F ) - 6(t' - t) 6(e F - e k )] , (21) 



corresponding (in frequency space) to simple poles, just like a Hartree Green's function. It 
is well known that the Kohn-Sham single-particle eigenvalues e k are not physical except at 
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FIG. 2: Equation for the density, showing the equivalence of the full and Kohn-Sham densities. 





FIG. 3: Cancellation of the density-density correlator with 5 Jo/ 5 p. 



the Fermi surface 0, Nevertheless, the trace of this Green's function gives the complete 
ground-state density p(x) (that is, the exact result if we calculate to all orders). 

We can easily show diagrammatically that Eq. ()18j) implies that the density obtained from 
the Kohn-Sham Green's function is, as advertised, exactly equal to that obtained from the 
exact Green's function. The density can be directly expressed in terms of the Kohn-Sham 
Green's function with equal arguments as 



p(x) = -iG 



ks 



— iuG^ s {x. .r 



(22) 



where v is the spin-isospin degeneracy. In Fig. ^ we have rewritten the last term in the 
Dyson equation (fTH|) for the exact Green's function using 



5T; 



5T int 5J (z) 

8p(y) J $M Z ) s p(y) 



d 3 z 



(23) 



where 8Jo(z)/5p(y) = [5 2 Wo/5Jo(y)5Jo(z)]~ 1 , which is minus the inverse density-density 
correlator Ilq. is represented with a double line (with no arrow). The result of carrying 
out Eq. on Eq. (jTSjl is shown in Fig. |21 where the last two diagrams cancel as in Fig. 01 
Note that while similar cancellations were shown in Ref. in the special case of zero-range 
interactions, the result here is completely general. Thus we see that the exact density is 
reproduced by the Kohn-Sham Green's function by construction. 

To illustrate some issues in comparing Kohn-Sham and exact Green's functions, we apply 
the formalism with the effective Lagrangian for dilute Fermi systems used in prior investi- 
gations: 



id t + 



2MJ 



16 



a 



(24) 



where V = is the Galilean invariant derivative and h.c. denotes the Hermitian 

conjugate. To describe trapped fermions, we add to the Lagrangian a term for an external 
confining potential v (x) coupled to the density operator v(x)if)ty [l4T |. For the numerical 
calculations, we take the potential to be an isotropic harmonic confining potential, 



v x 



2 I |2 

-moo x 



(25) 
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although the discussion holds for a general non-vanishing external potential. 

We repeat the previous development to introduce a second energy functional with an 
additional local source coupled to the kinetic energy density, following Ref. j 1 31 ] - The com- 
parison of results from the two functionals illustrates how the Kohn-Sham single-particle 
spectrum can be significantly different even though the bulk observables are essentially equal 
la l. So we consider 



and the corresponding effective action 

r>, r, £] = W[J, v,t]-J d A v J(y)p(y) - f d 4 y vivMv) , (27) 

with kinetic energy density 

t(x) ee (Vtf(x) ■ ViKx)) Jttl> t = &W }^ ] ■ (28) 

(We use superscript primes on the functionals, and ks' on the self-energies and Green's 
functions to distinguish the following quantities from those without rj or r dependence.) 
Each step goes through with straightforward generalizations, yielding Eq. (J18)) again, but 
now with 

yU = ^nJp,^] + + v ^ m v ^ p _ 

= E' w (y', VU + [My') + • V yVo (y')} 6 yS S 4 {y' - y) (29) 

after J(y') and r)(y') are set to zero. [Note that the gradients act on the G ks ^s to produce r 
after partial integrations in Eq. 

These two functionals were compared in Ref. 0] for a dilute gas of fermions in a harmonic 
trap. Two sets of parameters were used to illustrate the impact of a larger effective mass 
M*(x), which appears only in the "pr" (primed) formalism. Even though the fermion density 
and energy per particle for the p and pr functionals were very similar, the single-particle 
spectra have significant and systematic differences (see Ref. jl3( for details and figures). 
We can understand the systematics of the difference by comparing Kohn-Sham and exact 
spectra for a uniform system. We will drop the non-Hartree-Fock terms, which have been 
treated in LDA in both cases and which contribute equally to the energy spectra. We note 
that the terms in the pr functional correspond directly with terms in conventional Skyrme 
energy functionals [Kl- 
in the p case, the Kohn-Sham equation for the single-particle orbital [l^ (with external 
potential set to zero) leads to the spectrum 

= (30) 



where 
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In the pr case, we find a different spectrum 



where 



and 



k 2 



v-1 „ (v- 1 C 2 , v+ 1 C 2 



J °=~— c °»~ — f + — T>'> (33) 



1 1 (v-\C 2 v + \C' 2 



Vo = WYJ+ —-F + — -T P- (34) 



2M* 2M ' 2M V i/ 4 i/ 4 
Using r = |fcpp, the difference in the spectra simplifies to 

-P _ ( v — 1 C 2 V + 1 CgX 2 ,2\ 



Thus, the spectra differ for all momentum states except at the Fermi surface, where the 
spectra coincide as expected in Kohn-Sham DFT. In detail, the pr spectrum includes explicit 
momentum dependence that is converted to density dependence (i.e., kp dependence) in the 
p spectrum. We can also compare the Kohn-Sham spectra to that of the Green's function 
in the Hartree-Fock approximation, where we find that the e pT spectrum is the same as the 
Hartree-Fock spectrum. Indeed, for this approximation the Jo and 7/o terms in Eq. (|29|) 
precisely cancel against S^,. In contrast, Eq. (fT^j) yields a net contribution that shifts the 
Kohn-Sham spectrum to the Hartree-Fock spectrum. 

This example illustrates how individual Hartree-Fock self-energies in a gradient expan- 
sion can be completely included by adding the corresponding source terms. (A different 
example with covariant energy functionals is given in Ref. jll|.) The exact cancellations are 
only possible for local self-energies, which means Hartree-Fock. Beyond Hartree-Fock, the 
single-particle spectrum from the Kohn-Sham and exact Green's functions will necessarily 
differ. We can anticipate that self-energies with large non-localities will lead to the most 
significant differences. This is consistent with the expectation that low-lying vibrational 
states can account for the difference in level density between Skyrme (or other mean-field) 
and experimental spectra near the Fermi surface [2J, |5( . 

In this work, we have illustrated the relationship between Kohn-Sham and exact Green's 
functions within an effective action formalism. This approach goes beyond the observation 
that single-particle properties are not reliably calculated in terms of Kohn-Sham orbitals 
and eigenvalues. The formalism presents two ways to improve single-particle spectra. The 
Kohn-Sham spectra became closer to the exact spectra with the addition of appropriate 
sources. It is tempting to conclude that adding additional sources can always improve the 
Kohn-Sham single-particle spectrum, but this will require tests beyond the Hartree-Fock 
level. More generally, Eq. (|T5j) shows how to calculate single-particle quantities in terms 
of Kohn-Sham propagators at the same level of approximation (which is determined by 
the truncation of T\ nt ). In future work, the formalism will be applied to the calculation 
of spectral functions and the effect of low-lying vibrational states on the spectra tested by 
including self-energy diagrams that sum particle-hole bubbles. 
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